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Abstract 



We study the stability of the reconstruction of the scattering and absorption 
i^ ■ coefficients in a stationary linear transport equation from knowledge of the full 

albedo operator in dimension n > 3. The albedo operator is defined as the map- 
"^ I ping from the incoming boundary conditions to the outgoing transport solution 

J2 ' at the boundary of a compact and convex domain. The uniqueness of the recon- 

struction was proved in [H |4] and partial stability estimates were obtained in [13] 
for spatially independent scattering coefficients. We generalize these results and 

prove an L^-stability estimate for spatially dependent scattering coefficients. 
> 

o 

m ■ 1 Introduction 

^ I Let the spatial domain X C M", tt, > 2, be a convex, open bounded subset with C^ 

boundary dX, and let the velocity domain V be S"~^ or an open subset of M" which 

O ■ satisfies inf„gy \v\ > 0. Let T± = {{x, v) G dX x V; ±n{x)v > 0} where n{x) denotes the 

outward normal vector to dX at x G dX. The set r_ is the set of incoming boundary 

k> ! condition while r_|_ is the set where we measure the outgoing solution to the following 

H ' stationary linear Boltzmann transport equation in X x V^: 



vVxf{x,v) + a{x,v)f{x,v)— / k{x,v' ,v)f{x,v')dv' = in X x V, (LI) 

Jv 

/|r_=/- 

Here, f{x,v) models the density of particles at position x E X with velocity v eV. 
The albedo operator A is then defined by 

A:f-^f\r^, (L2) 
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where f{x,v) is the solution to (II. ip . The inverse transport problem consists of re- 
constructing the absorption coefficient a{x,v) and the scattering coefficient k{x,v',v) 
from knowledge of A. Stability estimates aim at controlling the variations in the recon- 
structed coefficients a{x, v) and k{x, v', v) from variations in A in suitable metrics. 

The forward transport equation has been analyzed in e.g. [5], [6l [8]. The inverse 
transport problem has been addressed in e.g. [21 HI [9l [10] with stability estimates 
obtained in [9l [13]. For the two-dimensional case, in which proofs of uniqueness of the 
scattering coefficient are available only when it is sufficiently small or independent of 
the spatial variable, we refer the reader to e.g. [H [HI [12]. 

To obtain our stability estimates, we follow a methodology based on the decompo- 
sition of the albedo operator into singular components [3l H] and the use of appropriate 
functions on r± with decreasing support [13]. In dimensions n > 3 the contribution 
due to single scattering is more singular than the contribution due to higher orders of 
scattering. As a consequence, the single scattering in a direction v' generated by a delta 
function /_ = 6xg{x)S{v — Vq) is a one- dimensional curve on dX. In order to obtain 
general stability estimates for the scattering coefficient, one way to proceed is to con- 
struct test functions whose support converges to that specific curve. It turns out that 
it is simpler to work in a geometry in which this curve becomes a straight line. 

We now briefly introduce that geometry and refer the reader to section [2] below for 
a formal presentation. Let i? be a positive real constant such that X is included in 
the ball B{R) of radius R centered at x = 0. On {B{R)\X) x V, the absorption and 
scattering coefficients vanish and we may solve the equation vVxf = 0. This allows us 
to map back the incoming conditions /„ on r_ as incoming conditions, which we shall 
still denote by /_, on F_ and map forward the outgoing solution /|r+ to an outgoing 
solution /+ on F^, where we have defined 

F± := {(x ± Rv, v) eW" xV for {x, v) eW xV s.t. vx = 0, \x\ < R}. (1.3) 

In other words, F± is the union for each v E V oi the spatial points on a disc of radius 
i? in a plane orthogonal to v and tangent to the sphere of radius R. 

The incoming boundary condition is thus now defined on F_ while measurements 
occur on F+ and we may define the albedo operator still called A as an operator mapping 
/_ defined on F_ to the outgoing solution /+ on F+. We may now verify that the single 
scattering in a direction v' generated by a delta function /_ = S{x — Xo)6{v — Vq) for 
(xo, Vq) G F_ is a one-dimensional segment in F+; see Fig. [H Note also that the geometry 
we consider here may be more practical than the geometry based on r±. Indeed, we 
assume that the incoming conditions are generated on a plane for each direction of 
incidence, and, more importantly, that our measurements are acquired on a plane for 
each outgoing direction. This is how the collimators used in Computerized Tomography 
[7] are currently set up. 

Under appropriate assumptions on the coefficients, we aim to show that ^ is a 
well posed operator from L^{FJ) to L^(F+). We shall then obtain a stability esti- 




Figure 1: Geometry of the ballistic and single scattering components in dimension n = 3. 
The source term is non-zero in the vicinity (in F_) of Xq — Rv'^ in Pi = —Rv'q + H^/ {R). 
The ballistic part is non-zero in the vicinity (in F+) of x'q + Rv'q in P2 = ^^0 + ^v' (R)- 
The thick "line" represents the support of the single scattering contribution in the 
vicinity (in F+) of the segment {x'q — sVq + t{s)v; s G (— -R, R)} C P3 = Rv + Ily{R). See 
text for the notation. 



mate for the reconstruction a{x) (or a{x, \v\)) and k{x,v,v') with respect to the norm 
£(Li(F„),Li(F+))ofA 

The rest of the paper is structured as follows. Because our geometry is not stan- 
dard, we present a detailed analysis of the linear transport equation and of the singular 
decomposition of the albedo operator in section [21 Most of the material in that section 
is similar to that in [3]. One of the main physical constraints in the existence of solu- 
tions to (II. ip is that the system be "sub critical", in the sense that the "production" of 
particles by the scattering term involving the scattering coefficient k{x,v,v') has to be 
compensated by the absorption of particles and the leakage of particles at the domain's 
boundary. Although this may be seen implicitly in [3], we state explicitly that the de- 
composition of the albedo operator used in the stability estimates holds as soon as the 
forward transport problem is well-posed in a reasonable way. 

The stability results are stated in section [31 Under additional continuity assumptions 
on the absorption and scattering coefficients, we obtain that (i) the exponential of line 
integrals of the absorption coefficient and (ii) the scattering coefficient multiplied by the 
exponential of the integral of the absorption coefficient on a broken line are both stably 
determined by A in £(L^(F_), L^(F+)); see Theorem 13.21 Under additional regularity 
hypotheses on the absorption coefficient, we obtain a stability result for the absorption 



coefficient in some Sobolev space H'^ and for the scattering coefficient in the L^ norm. 
The stabihty results in the geometry of (ll.ip are presented in section HI The proof of the 
stabihty results and the construction of the appropriate test functions are presented in 
section [5l Several proofs on the decomposition of the albedo operator and the uniqueness 
of the transport equation have been postponed to sections [6] and [3, respectively. 

2 Transport equation and albedo operator 

We now state our main results on the stationary linear transport equation and the 
corresponding albedo operator. 

Let i? be a positive real constant and let n G N, n > 2. Let V be S"~^ or an open 
subset of M" which satisfies vq := inf^gy 1^;! > 0. For v E V , we define v := -rr. Then, 
we consider the open subset O oiMJ^ xV defined by 

O := {(x, i)) e M" X 1/ I \xv\ < R, \x- {xv)v\ < R}, (2.1) 

and let F be the set 

F:={{x,v)eWxV \xe n^(i?)}, (2.2) 

where 

Uy{R) := {x eW \ xv = 0, \x\ < i?}, (2.3) 

for aA\ V E V. For all v E V we also consider 

Rv + n^{R) ■.= {Rv + x\x eli^iR)}. (2.4) 

When V = S""'^, then F is an open subset of T§"~^ := {{x,v) G R'^xS'^-^ | vx = 0}, 
the tangent space to the unit sphere. When V is an open subset of M"" (which satisfies 
Vq = inf^gy |f I > 0) then F is an open subset of the 2n — 1 dimensional manifold 
{{x,v) eR"- X V \ vx = 0}. We also define F± by 

F± ■={{x±Rv,v) gM" X V\{x,v) gF}, (2.5) 

and recall that F_ is the set of incoming conditions for the transport equation while F^ 
is the set in which measurements are performed. 
We consider the space L^{0) with the usual norm 



o 



f \f{x,v)\dxdv, for / G L\0). (2.6) 

Jo 



We also consider the space L^ (F) defined as the completed Banach space of the vector 
space of compactly supported continuous functions on F for the norm 

II/I|f:= / / \v\\f{x,v)\dxdv, feL\F), (2.7) 



and similarly the spaces L^{F±) defined as the completed Banach space of the vector 
space of compactly supported continuous functions on F± for the norm 

||/±||f±:= / / \v\\U{x±Rv,v)\dxdv,UeL\F±). (2.8) 

Jv Jn^{R) 

We assume that: 

< a e L°°(M" X V), 

k{x, v' , v) is a measurable function on M" x 1/ x y, 

a(x, v) = k{x, v', v) = for (x, v', v) eW xV xV, \x\ > R, (2.9) 

< k{x, v', .) e L\V) for a.e. (x, v') eW' xV and 

ap{x,v') = jyk{x,v',v)dv belongs to L°°(M" x V). 

Under these conditions, we consider the stationary linear Boltzmann transport equation 

vV^f{x,v) + a{x,v)f{x,v)- k{x,v',v)f{x,v')dv' = in O, (2.10) 

Jv 

f\F. = f- 

Throughout the paper, for ?n G N and for any subset U of M™ we denote by xu the 
characteristic function defined by ^^/(a;) = if a; ^ f/ and Xu{-c) = 1 ii x E U. 

We now analyze the well-posedness of (12.101) . The following change of variables is 
useful. 

Lemma 2.1. For f e L^{0), we have: 

f{x,v)dxdv= j j / f{y±tv,v)dtdydv. (2-11) 

O Jv Jliv{R) J-R 



Proof. First using (12. ip we have 



f{x,v)dxdv= / / xo{.X:v)f{.x^v)dxdv. 
o Jv Jr" 

Then for a.e. v E V using the change of variables {x E M" | \v.x\ < R, \x — {vx)v\ < 
R} -^ Ilj,{R) X {—R, R), X i-^ {x — {vx)v, ±vx), we obtain (12. lip . D 

We introduce the following notation: 

To/ = —vVxI ill the distributional sense, Aif = —af, (2.12) 

A^f = f k{x, v', v)f{x, v')dv', Ti = To + Ai, T = To + Ai + ^2 = Ti + ^2, 
Jv 

and the Banach spaces 

W := {fEL\0,\v\dxdvy,TofEL^iO)}, ||/||w = llTo/Ho + |||t^|/||o, 
W := {fEL\Oy,TofEL\0)}, ||/|U = ||To/||o + ||/||o. 



We consider the space L{F±) defined as the completed Banach space of the vector space 
of compactly supported continuous functions on F± for the norm 

II/±I|l(f±) ■■= [ [ \f±{x'±Rv,v)\dx'dv, for /± e L(F±). (2.13) 

Note that W C W and L^{F±) C L{F±). The spaces VV and L{F±) are used only to 
define the unbounded operators T and Ti below. We obtain the following trace result. 

Lemma 2.2. We have 

\\flFjF,<C\\f\U (2.14) 

for f G W, where C = max((2i?)~-'^, 1) and 

||/|fJ|l(f±) < C'll/IU, (2.15) 

for few, where C = max((2i?)~\ v^^). 

Proof. Let / be a C^ function in M" x y with compact support. Then from (12.111) . it 
follows that 

ll/llw= / / / {\v\df{x' ±tv,v)\ + \v\\fix' ±tv,v)\)dtdx'dv. (2.16) 

Jv Jn^(R) J-R "t^ 

Let V E V and x' G M" such that vx = 0. Note that f{x' =F Rv,v) = f{x' ± ^-0,^) — 
Ilnifi^' ± sv,v)ds for all t G {-R,R). Hence \f{x' T Rv,v)\ < \f{x' ± tv,v)\ + 
f-R li^fi^' ^ ■^'^' v)\ds. Upon integrating the latter equality, we obtain 

\f{x'TRv,v)\<^l^ \f{x'±tv,v)\dt + j^ \j-f{x'±sv,v)\ds. (2.17) 

Combining (KW\i et (12X7]) . we obtain (1214]) . The proof of (|2:T5D is similar. D 

For a continuous function /_ on F_, we define the following extension of /_ in O: 
Jf^{x,v) = e-^^^'' ^^^^^' ''^^-'^'"'>'^' f^{x - {xv + R)v,v), {x,v) G O. (2.18) 

Lemma 2.3. For /_ G ^H^-) ^f'^^^ C* = 2/2(1 + Vfj-^HcrHoo), we /ia^e.- 

\\Jf-\\w<C\\f4F_. (2.19) 

Proof. Let /_ be a compactly supported continous function on F_. From (12. lip and 
fl2:T8D it follows that 

|||^|J/„||o=/ [\v\\f4x' -Rv,v)\ /e-l"l"'^-«"("'+^'''")^^rftrfx'(it;<2i?||/_||^_. (2.20) 
JvJn-u{R) J-R 

One can check that J/_ satisfies To J/_ = —AiJf^ in the distributional sense. Therefore 
using also (^M> we obtain ||ToJ/_||o + |||f|J/_||o < (1 + ll^ikTi) x lll^l-^/-l|o = 
2i?(l + V(^"'^||cr||oo)||/-||F_5 which proves the lemma. D 
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2.1 Existence theory for the albedo operator 

We consider the following unbounded operators: 

Ti/ = TJ, T/ = Tf, D(Ti) = D{T) = {/ G W ; f\F_ = 0}. (2.21) 

The operator Ti : -D(Ti) -^ L^{0) is close, one-to-one, onto, and its inverse Ti^^ is 
given for all / G L}{0) by 

fR+xv 

T^-'f{x,v) = -\v\-^ e-l"l"'^o-(--«'-)'^7(x-t{),^;)(it, {x,v)eO. (2.22) 

Jo 

Lemma 2.4. The following statements hold: 

i. The bounded operator \v\Ti^^ in L^{0) has norm less or equal to 2R and the 
bounded operator A2\v\~^ in L^{0) has norm less than \\\v\~^ap{x,v)\\L°^(o)- 

a. Under the hypothesis 

a - ffp > 0, (2.23) 

the bounded operator ^2X1"^ in L}{0) has norm less than 1 — e~^^^o IkplU. 

Hi. Assume either condition ( 12.23^ or 

2R\\\v\-^ap{x,v)\\L^^o) <1- (2.24) 

Then I + ^2X1" is invertible in L^{0). 

Lemma 12.41 is proved in section 7. We denote by K the bounded operator in 
L\0, \v\dxdv) defined hj K = Tr^^: 

Kf{x,v) = -\v\~^ / 6-^"^'' ^o-(^-'^'^)'i'(A2f)ix - tv,v)dt, {x,v) G X X V, 

Jo 

for all / G L^{0, \v\dxdv). The operator K also defines a bounded operator in L^{0). 
This allows us to recast the stationary linear Boltzmann transport equation as the 
following integral equation: 

{I + K)f = Jf^. (2.25) 

The existence theory for the above integral equation is addressed in the following result. 

Proposition 2.5. The following statements hold: 

i. The conditions (I2.26P and (I2.27P below are equivalent. 

The bounded operator I + K in L^{0) admits a bounded 
inverse in L^{0). (2.26) 

The bounded operator I + ^42X1""^ in L^{0) admits a bounded 
inverse in L^{0). (2.27) 



ii. Assume either (12.231) or (12.241) . Then condition (12.261) is satisfied. 
Hi. If (I2.26P is satisfied then 

the bounded operator I + K in L^{0, \v\dxdv) 

admits a bounded inverse in L^{0, \v\dxdv). (2.28) 

Proposition 12.51 is proved in Section 7. The following proposition deals with the 
existence of the albedo operator. 

Proposition 2.6. Assume (12.281) . Then 

i. the integral equation (I2.25P is uniquely solvable for all /_ G V-{FJ) and f eW ; 

ii. the albedo operator A : f- ^^ f+ = f\F^ is a bounded operator A : L^{F^) -^ 
Proposition 12.61 is proved in Section 7. 

2.2 Singular decomposition of the albedo operator 

We assume that condition (12.281) is satisfied. Let us consider the operator 
7^ : L\0, \v\dxdv) -^ L\F+), defined by 

^^7^^:=(J^2^)|^^, (2.29) 

for ip G L}{0, \v\dxdv). Using the equality (in the distributional sense) T^Kf = —Aif — 
A2/ for / G L^{0, \v\dxdv) and the boundedness of the operators Ai and A2 from 
L^(0, \v\dxdv) to L^{0, dxdv) and using (I2.14p . we obtain that 7^ is a well defined and 
bounded operator from L^{0, \v\dxdv) to V-{F^). We shall use the following lemma for 
the kernel distribution of TZ. 

Lemma 2.7. We have the following decomposition: 

n^l){x,v)= f3{x,v,x',v')ip{x',v')dx'dv', (2.30) 

Jo 

for a.e. {x,v) G F+ and for any ip G L^iO, \v\dxdv), where 

< \v'\-^p G L'^iO, L\F+)). (2.31) 

In addition if k E L°°(M" x V x V), then for any e' > 0, 5 > 0, and any 1 < p < 
1 + ^^ there exists some nonnegative constant C{e',6,p) such that 




(j){x, v)P{x, V, x', v')\v\dxdv 

VJRv+U^R) L^{0,,y) 

<C{e',6,p)( f f \^{x,v)fdxdvY +e'mL^^F^), (2.32) 

\Jv JRv+n^{R) J 
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for any continuous compactly supported function 4> on F+ such that suppi/) C {{x^v) G 
F+ I \v\ < S^^}, and where p'~ + p~'^ = 1. 

Lemma 12.71 is proved in Section 6. The last inequality shows that the kernel of 
the second scattering operator 7?. is more regular than is indicated in ( 12.3ip . When 
V is bounded, then we can choose e' = in (12.321) . in which case we obtain that 
\v'\~'^P G L°°{0,LP{F+)) for 1 < p < ;^. This regularity is sufficient (while that 
described in (I2.3ip is not) to show that multiple scattering contributions do not interfere 
with our stability estimates. Taking account of Lemma 12.71 we have the following 
decomposition for the albedo operator. 

Lemma 2.8. Under condition (I2.28p . the following equality in the distributional sense 
is valid 

A(l)-{x,v) = a{x,v,x',v')(j)-{x' — Rv',v')dx'dv' 

JvJuAR) (2.33) 

+ / f3{x,v,x\v'){{I + K)-^J(t).){x' ,v')dx'dv' , 
Jo 

for a.e. {x,v) G F+ and for any C^ compactly supported function 0_ on F_, where 



a{x,v,x',v') = ai{x,v,x' ,v') + a2{x,v,x' ,v'), (2.34) 

-(xv')v 



a,{x,v,x',v') = e-l^'l"'/o'^^(--.-)'^^5„(^;')5^^ , ,(:r'), (2.35) 



/ „,/\ l„.|— 1 / ~\v\^^La{x—sv,v)ds-~\v'\^^Jg ^^ " o-{x—tv—sv',v')ds 



a2[x, v,x ,v ) = \v\ / e 

'0 

xk{x - tv,v\v)d^_^^_^^^_^-yyix')dt, (2.36) 

for a.e. {x,v) G F^ and {x',v') G F, and where (3 is given by (12.301) . 

Lemma 12.81 is proved in Section 6. The above decomposition is similar to that 
obtained in [3l H] except that the multiple scattering contribution is written here in 
terms of the distribution kernel of TZ rather than that of TZ{I + K)~^J. 

3 Stability estimates 

In this section, we give stability estimates for the reconstruction of the absorption and 
scattering coefficient from the albedo operator following the approach in [T3] . 

We assume that conditions (12. 9p and (12.280 are satisfied and that there exists a 
convex open subset X of M" with C^ boundary dX such that X C -8(0, R) := {x G 

M" I |x| < R} and 

the function < (7\xxv is continous and bounded in X x V, 

the function < k\xxVxv is continous and bounded in X x 1/ x F, (3.1) 

a{x, v) = k{x, V, v') = for X ^ X, {v, v') e V x V. 

9 



Let (5", A;) be a pair of absorption and scattering coefficients that also satisfy (12. 9p . 
(12.28p . and (I3.ip . We denote by a superscript ~ any object (such as the albedo operator 
A or the distribution kernels cij, i = 1,2) associated to (5", k). 

Let {x'q, v'q) G F such that the intersection oi X and the straight line {tv^+XQ \ t G M} 
is not empty. The point {xq — Rvq, Vq) G F„ models the incoming condition and is fixed 
in the analysis that follows. For e > let /^ G C^(F„) such that ||/e||F_ = 1, /e > 
and supp/g C {{x' — Rv',v') G F_ \ \v' — Vq\ + \x' — Xq\ < e}. Hence \v'\f^ is a smooth 
approximation of the delta function on F_ at [xq — Rv'q, v'q) as e ^ 0^ and is thus an 
admissible incoming condition in V'{FJ). The support of /^ is represented in Fig. [H 

For a.e. (x, ti) G F, t G M and v' E V let E{x, t, v, v') > be defined by 



E(x t V v') := e~'^' ^ J_j^<7(x-sv,v)ds-\v'\ i/^ '•'" "•"' a{x-tv-sv',v')ds ^ /•2^2) 

Replacing a by a in (13. 2p we also define E(x,t,v,v') for a.e. {x,v) E F, t E M. and 
v' eV. 

Let (5 > and let (j) be any compactly supported continuous function on F+ such 
that \\4>\\oo < 1 and 

SUPP0 C {{x,v) E F+ I \v\ < S'^}. (3.3) 

Then using (I2.33P and (13.21) we obtain for 5 > that 




\v\(f){x,v){A- A)fs{x,v)dxdv = Ii{(p,e) + J2(0, e) + /3(0,e), (3.4) 



where 



I^{(j),e) = \v\(P{x + Rv,v)(e-\''\~'^^n-(---'-)'i^ 

Jv Jn^{R) ^ 

_^-\v\-' !^^a[x-sv^v)ds\ j^^^ _ fly^^^dxdv (3.5) 

12(0,5) = / / (f){x + Rv,v) / {k{x — tv,v',v)E{x,t,v,v') 
Jvxv Jn^{R) J~R 

—k{x — tv,v',v)E{x,t,v,v')] (3.6) 

/e(x — t-u — (x — tv)v' — Rv', v')dtdxdvdv' , 
h{<P,e) = 4(0,£)-/|(0,£), (3.7) 



and where 



Il{(f),e)= / / \v\(p{x,v) / (3{x,v,x',v'){{I + K)-^Jf,){x',v')dx'dv'dxdv, (3.8) 
iy Jrv+UvIr) Jo 

Il{(t),e)= f [ \v\(j){x,v) f /3{x,v,x',v'){{I + k)~^Jfe){x',v')dx'dv'dxdv. (3.9) 

Jv JRv+ny{R) Jo 
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In addition using the estimate 
of fe, we obtain 



< 1, item ii of Proposition 12.61 and the definition 




|v|(^(a;, v){A — A)fs{x, v)dxdv 

V JRv+Uy{R) 
< \\A- A\\c{L^F^),L^F+))\\fe\\F^ 



<\\iA-A).MW, 

•^- ^\\c{L^F^),L^F+))■ (3.10) 



3.1 First stability estimate 



We now prove a stabihty estimate under conditions (12. 9p . (12.281) . and (13. ip . Taking 
(I3.5p - (I3.7I) into account, we obtain the following preparatory lemma: 



Lemma 3.1. Assume that {(J,k) and (a, /c) satisfy conditions (12. 9p . (I2.28p . and (13. ip . 
Then the following limits and estimate hold: 



/i(0, e) — ^ 0(xo + Rv'f^, v'q) 



X e 



kol V-flO"(^o-H''"o)'^'^ _ g-l'^ol ^ f-R^i^'o-^'"'oA)'^^ 



,e) -^ 4(0) + J2^(0), 



for any compactly supported continuous function (p on F^, where 



/2(</') 



1 




{k-k){xQ + t'v'Q,VQ,v) 



JQ] Jv J-R 



X [(p{x + Rv, v)E{x, t, V, v'q)] t^t{=o'^y^,t',v) dt'dv, 



X^X{Xr^,Vf. ,t ,v) 



(3.11) 



(3.12) 



(3.13) 



R 



Iii<P) = VT] k{x', + t'v'„v'„v) 

[■yol Jv J-R 



X 



[x + Rv, v) {E - E) (x, t, V, v'o) ,^,w ,.,,,) dt'dv, 



t{x'^,v'fj,t' ,v) 



(3.14) 



where E and E are defined by (13. 2p and 

{t{x'o, v'q, t', v),x{x'o, v'o, t', v)) = (-{x'q + t'v'Q)v, x'o + t'v'o - {{x'o + t'v'Q)v)v ) , (3.15) 



for t' G M. In addition, for all e' > 0, 5 > and for all 1 < p < 1 + ^^ there exists 
some nonnegative real valued constant C{e',6,p) such that 

sup |/3(0,£:)| < C C(£',5,p) ( / / Xsnpp<pix + Rv,v)dxdv\ + e' ] , (3.16) 
£>o \ \Jv JuviR) / 
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for any compactly supported continuous function on F^, which satisfies \\4>\\od < 1 d'nd 
(13 -Sp for S > 0, where p'~ + p~^ = 1 and 

C := 2R\\{I + Ky \\c{L^o,\v\dxdv)) + '^R\\iI + Ky \\c{L^io,\v\dxdv))- (3-17) 

Lemma [Hi] is proved in Section 5. 

Taking account of Lemma [3.11 and (13.101) . and choosing an appropriate sequence of 
functions "</>" , we obtain the main result of the paper: 



Theorem 3.2. Assume that n > 3 and that (cr^k) and (a, /c) satisfy conditions fl2.9p . 
(12.281) . and (13. ip . Then the following estimates are valid: 



K\ ^ I-R(^i^'o-'^'"'o''"'o)'^'^ 



K\ V-flO"K-'^^o'''o)"'^ 



< \\A-A\\cimF.),L 



H^+))5 



(3.18) 



.' 1-1 




R 



'V J-R 

J l-lll;:; 



[k - k){xQ + t'vQ, v'q, v) [E{x, t, V, Vq)] t=t(,^,„^,t',„) dt'dv 



< 2R\v'q\ ||^p(a;o + t''yo'^o)IU^(iR,,) sup {E - E){x,t,v,v'f^) 



(x,v)eF 

tgi 



\A-A\ 



C{LHF-),L^{F+)), 



(3.19) 



where E and E are defined by (13. 2p . and where {t{xQ,v'f^,t',v),x{xQ,v'Q,t',v)) is defined 
by ([3J[5P for t' eR andv e V. 

Theorem 13.21 is proved in Section 5. 



Remark 3.3. Under condition (I2.26p . we can obtain similar estimates to those in The- 
orem \3^ for the albedo operator defined on L{F) with values in L{F+). Note that 
L{F±) = L^{F±) when V is bounded. 



3.2 Stability results under additional regularity assumptions 

The second inequahty in Theorem 13.21 provides an L^ stabihty result for k{x,v',v) pro- 
vided that a{x,v) is known. The first inequality in Theorem 13.21 shows that the Radon 
transform of a{x, v) is stably determined by the albedo operator. Because the inverse 
Radon transform is an unbounded operation, additional constraints, including regular- 
ity constraints, on a are necessary to obtain a stable reconstruction. We assume here 
that 

{^ \veV] = §"~\ Vo := sup \v\ < oo, (3.20) 

\V\ v€V 
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and that the absorption coefficient a does not depend on the velocity variable, i.e. 
a{x^v) = cr(x), X G M""; see also remark [331 below. Then let 

M := {{a{x),k{x,v',v)) G L°°(M") x L°°(M" x V x V) \ {a,k) satisfies ([31]), 
(12:91) and ^M^, and a|x G //^+'(X), \\a\\^^+,^^^ < M, \\ap\\^ < M}, (3.21) 

for some f > and M > 0. Using Theorem 13.21 for any (xQ,i'o) € -F such that the 
intersection of X and the straight line {tv^ + Xg | t G M} is not empty, we obtain the 
following theorem. 

Theorem 3.4. Assume that n > 3. Under condition (13.201) . for any {a,k) G Mi and 
[a, k) E M. the following stability estimates are valid: 

Ik - 5-||//.(x) < Cill^ - ^||£(Li(F_),Li(F+))' (3.22) 

where -\<s<l + f,e = ^^±^5?, and C^ = C,iR,X,Yo,Vo, M,s,r); 



2 ' ' n+l+2f 




V J-R 



k{xo + t'v'o, Vq, v) - k{xQ + t'v'o, v'q, v) 



dt'dv 



< C2M - '^\\c(L^F^),L^F+)) [^ + M - ^\\c{L^(F-),L'^{F+))) ' (3.23) 

for {xq,Vq) G F such that Xq + t'v'^ G X for some t' G M, and where 6 = J^^^lf ' 
< r < f, and C2 = C2{R, X, Vq, Vq, M, r, f); in addition, 

\\k — A;||Li(]8nx\/xy) 

where 9 = ^^j^, <r <f, and C3 = C3(/2,X, vq, Vq, M, r, f). 

Theorem 13.41 is proved in Section 5. 

Remark 3.5. Theorem \3.4\ can be extended to the case a = a{x, \v\) and V = {v E 
M" I < Ai < |f I < A2 < 00}. In this case the class M. is replaced by the class 

N := {{a{xM),k{x,v\v)) eL'^iW' x V) x L°°(M" x V xV)\ 

(ct, A;) satisfies ([31]), ([23]) and ^TM . ||ap||oo < M and for any A G (Ai, A2), 

a,x(.,A)Gift+^(X), sup ||a(.,A)||^n+. <M}. (3.25) 

AG(Ai,A2) 

Then the left-hand side of (13.221) is replaced by ^v^Y^x^iXiM) lk(-' '^)ll//f+*'fx) 

whereas the right-hand side of ( 13. 22^ and estimates ( 13. 23^ - (13.240 remain unchanged 



(see the proof of Theorem 3.4) 
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4 Stability in r± 

We now come back to the original geometry in (11.11) and present a similar stability result 
(Theorem 14.31 below) to Theorem 13. 4[ The case of a scattering coefficient k{x, v', v) = 
k{v',v) that does not depend of the space variable x was studied in [13]. We now 
introduce the notation we need to state our stability result. 

Recall that X C M", n > 2, is an open bounded subset with C^ boundary dX, and 
that V is S"""*^ or an open subset of M" which satisfies vq := inf„gy \v\ > 0, and that the 
linear stationary Boltzmann transport equation in X x y takes the form 



vVxf{x,v) + a{x,v)f{x,v)— / k{x,v',v)f{x,v')dv' = OmXxV, (4.1) 

Jv 
/|r_=/- 

We assume that {a, k) is admissible if 

0<aeL^{X xV), 

k{x, v' , v) is a measurable function on X x y x y, and ,. ^ 

< k{x, v\ .) G L\V) for a.e. (a;, v')eXxV ^ ^ 



ap[x,v 



Jyk{x,v',v)dv belongs to L°°{X x V). 



For (x, v) G {X xV)U T^, let t±{x, v) be the real number defined by r±(x, v) = sup{t > 
\x ± sv ^ X for all s G (0, t)}. For {x, v) ^ X x V, let r(x, v) be defined by 

t{x, v) = T^{x,v) + T^{x,v). 

For (x,v) G Fzp, we put t{x,v) = t±{x,v). We consider the measure dC,{x,v) = 
\n{x)v\dii{x)dv on F-|-. We still use the notation Tq, Ti, T, Ai, and A2 as in (12.121) 
and introduce the following Banach space 

W := {/ G L\X X V)-Tof G L\X x V), r'^ f G L\X x V)}, 

ll/llw = WTofWi^i^xxV) + Ik /llLi(Xxy)- 

We recall the following trace formula (see Theorem 2.1 of [1]) 

||/|rJ|Li(r±,dO < ll/llw, for/GW. (4.3) 

Estimate (14. 3 p is the analog of the estimate (I2.14p in the previous measurement setting. 
For a continuous function /„ on r_, we define J f- as the extension oi f_ m X xV 
given by : 

Jf^{x,v) = e~ fo'^^'"^ -i^-'^'^)d^ f^^x - T^{x,v)v,v), {x,v) eX xV. (4.4) 
Note that JT" has the following trace property (see Proposition 2.1 of [4]): 

||J/-||w<C||/„|Ui(r„,dO, (4.5) 

for /„ G L^{T^,d^), where C = 1 + diam(X)vg'"'^||a||oo and where diam(X) := 
sup^ygjf \x — y\. Estimate (14. 5 p is the analog of estimate (I2.19P in the previous mea- 
surement setting. 
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4.1 Existence theory for the albedo operator 

We denote by /C the bounded operator of L^{X x V, T~^dxdv) defined by 

/•T-{x,v) P 

}Cf{x, v) = - e- ^0 '^{=^-sv,v)ds / ^^^ ^/^ ^^^^^ _ ^^^ y'-^fiy'dt^ (x, v)eX xV. 

Jo Jv 

for all / G L^{X x V, T~^dxdv). We transform the stationary linear Boltzniann transport 
equation (14. ip into the following integral equation 

{I + IC)f = Jf^. (4.6) 

We have the following proposition, which is the analog of Proposition I2.6[ 

Proposition 4.1. Assume that 

the bounded operator I + K, in L^{X x V, T~^dxdv) admits a hounded inverse 
znL\X xV,T-^dxdv). (4.7) 

Then 

i. the integral equation (14. 6 P is uniquely solvable for all /_ G L^(r_, d^), and f G W, 

a. the operator A : L^(T_, d^) -^ L-'^(r+, d^), /_ -^ /|r+; is a bounded operator. This 
operator is called the albedo operator A for (14. ip . 

The above proposition can be proved by slightly modifying the proofs of Propositions 
2.3 and 2.4 of [1]. 

Remark 4.2. i. Assume that X is also convex. Let f G L^{F-i-) be such that 
supp/ C {{x,v) E F± \ X + tv E X for some t G M}^ where F± is defined by 
()2.5p and R > diam(X). Then we obtain that: 

/ / f{x ± Rv,v)\v\dxdv = / f{'j±{x,v),v)d^{x,v), (4.8) 

Jv in„(i?) Jr± 

where 7±(x, v) = x — {xv)v±Rv for any (x, v) eT±. Therefore, considering results 
on existence of the albedo operator A obtained in ^ and our assumptions (12. 9p . 
equality (14.80 leads us to define the albedo operator A from L^(F_) to L^(F+). 

a. The condition (14. 7p is satisfied under either of the following constraints: 

Ik^plloo < 1, (4.9) 

a-ap > 0. (4.10) 
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Hi. Assume that 

the bounded operator I + K, in L^{X x V) admits a hounded inverse 
mL\XxV). (4.11) 

Then we can define the albedo operator from L^(T^,d^) to L^(T^,d^) where d^ = 
min(T, X)d^ and where X is a positive constant. To prove this latter statement, we 
need trace results for the functions f eW := {f e L^{X xV) \ Tof E L^{X xV)}. 

iv. Under (14. 9 p and the condition ||ra|| < oo, the existence of the albedo operator 
A : L^(T_,d^) -^ L^{T+,d^) is proved in ^ (Proposition 2.3) when V is an open 
subset ofW^ (the condition inf^gy \v\ > zs not required). 

V. Under the condition a — ap > u > 0, the existence of the albedo operator A : 
L^(r_,d^) — > L^(T^,d^) is proved in ^ (Proposition 2.4) when V is an open 
subset o/M" (the condition inf^gy 1^;! > is not required). 

Finally under (14. 7p . we also obtain a decomposition of the albedo operator A similar 
to that of A given in Lemma 12.81 

4.2 Stability estimates 

We assume that X is convex and 

the function < a is continous and bounded on X x V, ,. . 

the function < A; is continous and bounded on X x V^ x V. 

Let (a, A;) be a pair of absorption and scattering coefficients that also satisfy (I4.12p . 
(14. 2 p and (14. 7p . Let A be the albedo operator from L^(r_, d^) to L^(r+, d^) associated to 
(5", A;). We can now obtain stability results similar to those in Lemma [3. II and Theorem 
13.21 Consider 

M := {{a{x),k{x,v\v))eH'^^^{X)xC{XxV xV)\ 

{a, k) satisfies glD, (iTD, M^^+,^^^ < M, \\ap\\^ < M} (4.13) 

for some f > and M > 0. We obtain the following theorem. 

Theorem 4.3. Assume n > 3. Under conditions (I3.20p . for any {a,k) G M and 
{a, k) G M, the following stability estimates are valid: 

\W - ^Wn^ix) < Ci\\A - A||^(^i(p_^^g)^^i(p_^^^^)), (4.14) 

yjhere -| < s < f + f, ^ = ^^gg^?, and C, = Ci(X, Vq, Vq, M, 5,f); 

k{xQ + t'v'f^,VQ,v) -k{x'f^ + t'vQ,v'f^,v) dt'dv (4.15) 




'V Jo 

< C2IIA - A||£(^i(p_^^^)_^i(r^_^^)) [I + II A - A||^^^i(p___^^^^^i(p_^^^^^-, 
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for (xq, v'q) G r_ and where 6 = ^1''-^ , 2^= , < r < r, and C2 = C2{X, vq, Vq, M, r, r). ^45 
a consequence, we have 

\\k - k\\L^(^xxVxV) 

< C'sllA - A||_c(^i(p_^^^)_^i(p^_^^)) (^1 + ||A - M\c{L^(r.,dO,LH^+40)) ' (4-16) 

where 9 = ^^y^, <r <f, and C3 = C3(X, vq, Vq, M, r, f). 
The proof of Theorem 14.31 is similar to that of Theorem 13.41 

Remark 4.4. Theorem \4-3\ can also he extended to the case a = a{x, \v\) and V = {v & 

W I < Ai < \v\ < A2 < 00}. In this case the class M is replaced by the class 

N := {{a{x,\v\),k{x,v',v))eC{X xV)xC{X xV xV) \ 

(cr, A;) satisfies (14.121) . (14. 7p . ||crp||oo < M, and for any A G (Ai, A2), 
a|x(.,A) G H^+'{X), sup ||a(.,A)||^n < M}. (4.17) 

Ae(Ai,A2) 

Then the left-hand side of (14.141) is replaced by s'^P\(z(^Xj^^x2) W'^i-^ ^)\\h^+^(x) 
whereas the right-hand side of (14.141) and estimates (14. 151) -( 14.161) remain 
unchanged. 

5 Proof of the stability results 

We now prove Lemma 3.1 and Theorems 3.1 and 3.2. 

Proof of Lemma \3.1\ Using the fact that X is a convex subset of M" with C^ boundary 
and using (13. ip . we obtain that 

the function 'RxMP xV 3 {t,x,v) ^>- I a{x — sv, v)ds is continous 

J-R 

at any point {i, x, v) such that x -\- r]v E X ioT some real 77. (5-1) 

The same statement holds by replacing a by a. From ( 13. 5p . it follows that 

Ii{(f),e)= / / ^i{x,v)fe{x - Rv,v)\v\dxdv, (5.2) 

•'^ •' \A<R 



where $1 is the bounded function on F defined for (x^v) G -F by 
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From (15.11) and the continuity of 0, it follows that $i is continuous at the point 
{x'q,v'q) in F. Therefore using (15.20 and the definition of the functions /e, we obtain 
lim£^o+ h{<P,£) = ^i{xq,Vq), which implies (13.111) . Performing the change of variables 
X — tv = x' + t'v' with x'v' = {"dtdx = dt'dx'") in formula (13.61) and using (13.11) . we 
obtain 



l2((f^,e) 




<^2,t',v{x'^^')fe{x' - Rv',v')\v'\dx'dv' | dvdt' , (5.4) 

RJV \JV JU^,{R) ' 



where 
and 



$2,t',.(a^',^^') = 0if x' + t'7;'^X, 



^2,t',v{x', v') = —- ( k{x' + t'v', v', v) [(f){x + Rv, v)E{x, t, v, v')] t=t(x' y ,t' ,v) 

\U \ \ x^x{x' ,v^ ,t^ ,v) 



(5.5) 



— k{x' + t'v', v' , v) (j){x + Rv, v)E{x, t, v, v') 
if x' + t'i' e X, 



t^t{x' ,v' ,t' ,v) 



(5.6) 



for (f, v) e {-R, R) X V, where 

(t(x', v', t', v),x{x', v , t', v)) := {-{x' + t'v')v, x + t'v' - {x' + t'v')v ) , (5.7) 



for x' E W, V, v' EV,t' E M. Let t' E {-R, R) such that x'q + t'v'^ E X, and let v eV. 
Then from (^^, ([SII])-(IS2D and (jSSD-dSZD it follows that ^2,t\v is continuous at the 
point {x'q,v'q). Hence 

^2,t',vix',v')fe{x' - Rv',v')\v'\dx'dv' ^ ^2,t',vix'Q,v'Q), as £ ^ 0+. (5.8) 




'V 



x'\<R 



Moreover using (I5.5p - (l5.6j) and using the estimate a > (and (13. 2p ) and the equality 
||/e||F_ = 1 we obtain 



'^2,t',vW^'"')feix' - Rv',v')\v'\dx'dv' 



x'\<R 



^ life + k\\oo\\(t>\\L^(F+)Xsnppy<t>{v) ,^ g^ 

~ Vo 



for {t',v) E {-R,R) X V, where snppyc/) = {v E V \ 3x E M", xi) = R, (j){x,v) ^ 0}. 
From (ESD, (ESI), (EH) it follows that \mi^^Q+ l2{(^,e) = j!!j^Jy^2,t',v'ix'Q,v'o)dv' dt', 
which implies ( 13.121) . 

It remains to prove ( 13.161) . We first estimate sup^^o 1-^3(05^)1- Using ( I2.28p . the 
estimate \\Jfe\\L^{o,\v\dxdv) < 2i?||/e||F_ and the equality ||/e||F_ = 1, we obtain 

II (/ + KT'JfeWmoMdxdv) < 2i?||(/ + Ky'\\ciLHo,\v\dxdv)), e > 0. (5.10) 
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Let e' > and 1 < p < 1 + ^ and p~^ + p'~^ = 1. Using (|23B . (ISlSD . (I2:32|) and 
the estimate ||^||l°°(f+) < 1 we obtain 



\lli<l^,e) 



{{I + K) ^ J fe){x\v') I (f){x,v)P{x,v,x',v')\v\dxdvdx'dv' 

o Jv Jrv+U^{R) 



<||(/ + /"i') Jf£\\L^{0,\v\dxdv) 



0(x, v)\v\——P{x, V, x' , v')dxdv 

V JRv+uJ{R) F I 




L°°(0) 



1 



<—Ci[C{e',6,p) I / Xsnpp^ix, v)dxdv] + e' 

vo \ \Jv JRv+U^{R) J 



where C(£',(5, p) is the constant from (12.321) . and 

C\ = 2R\\{I + Ry \\c{Li(o,\v\dxdv))- 



(5.11) 



(5.12) 



Replacing K, J, a and f3 by K, J, a and /3 in (I5.1ip - (l5.12p . we obtain an estimate for 
sup^^qI^{(J),£). Combining these estimates with (13.71) . we obtain (13.161) . D 

Proof of Theorem \3.^ Let ei > and let 0^^ be any compactly supported continuous 



£1 

function on F+ which satisfies < (^^^ < 1 and 

4>ei{x + Rv,v) = 1 for (x, v) G F, |a; — x'q\ + \v — Vq\ < — , 
suppc^ej C {(x, v) G F+, \x — Rv — Xq\ + \v — Vq\ < £i}. 

From fl3:TTD and (KT3\i it follows that 



lim lim Ii{e,(l)ei] 

£1^0+ £^0 + 



''ol ^ I-R'^i^'o-^^'o''"'o)d-'^ — p-Kl V-h'^(^0-H'*'o)'^'5 



From fl3A6D and flSAij) it follows that 



lim lim sup \l3{e,(j)ei)\ = 0. 

£1^0+ £^0+ 



(5.13) 
(5.14) 

(5.15) 
(5.16) 



From fl3:T3D . (137HD . it follows that 



2R. 



I^2(0£i) + ^2(0£i)l < (||A:||oo+ ||A:||oo) / Xsnppyc^,^{v)dv, 

Vo Jv 



(5.17) 



where suppy</)ej = {ti G F | 3x G M", xt) = R, (j)si{x,v) ^ 0}. Note that using (I5.14p . 
we obtain 



Xsupp^.0e, {v)dxdv < dv ^ 0, as £i ^ 0^ 

V J, "f|; 



(5.18) 
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From (I5A8D and (ISTTj) it follows that 



1^2(0..; 



+ /|( 



^ei, 



0, as £1 ^ 0" 



(5.19) 



Note also that from (13.41) and (I3.10p . it follows that 



|A(0£i,£)| < \\A- A\\c{L^^F.),L^{F+)) + \ll{4>e,^e)+ll 



^£11 



01 + 1/3(0.,, £)|, (5.20) 



for £ > and ei > 0. Combining (lOOD (with "0"= ^.J, (ISJ^) . (ISlTej) and (ISTTOj) . we 
obtain (13.18^ . This provides us with a stability result for the absorption coefficient. 

It remains to obtain a stability result for the scattering coefficient. We first construct 
an appropriate set of functions "0" (see (I5.28P below). The objective is to construct a 
sequence of such (smooth) functions whose support converges to the line in F+ where 
single scattering is restricted; see Fig{Tl Moreover, we want these functions to be good 
approximations of the sign of A; — A; on that support. This is the main new ingredient 
that allows us to obtain stability for spatially dependent scattering coefficients. More 
precisely, let U := {{t',v) e R x V \ x'q + t'v'Q E X and {k - A;)(x'o + tX,w^,w) > 0}. 
Using (13. ip . it follows that U is an open subset of M x F. Let (Km) a sequence of 
compact sets such that IJmeN ^m ~ ^ ^^^ ^m ^ -^m+i for ra E N. For m G N let 
Xm e C°°(M X V, M) such that xk^ <Xm<Xu, and let 



Pv 



^/Cm -L- 



(5.21) 



Thus using (13. ip we obtain 

lim {k-~k){x'o + t'v'f„v'o,v)p^{t',v) = \k-~k\{x'o + t'v'Q,v'o,v), (5.22) 

m— >+oo 

for V eV and t' E R such that x'q + t'v'^ E XU {W\X). For {x,v) E F+ such that v 
and v'q are linearly independent, we define 



d{x,v) :- 



X — x'q — ((x — Xq)v)v — ((x — x'q) 



V - (VaV)v'n , V - (vi,v)v'n 



JQVJUQ 



JQUJUQ 



1 - (v'qv) 



ln,\2 



I ,^\2 



WcV 



(5.23) 



For {x^v) E F+ such that v and v'q are linearly independent, we verify that d{x^v) 
infj^t/gK \x + tv — (x'q + ffp)! and the infimum is reached at 



(t, t') 



[x'q -x){v- {vv'q)v'q) {x - x[))(w^ - {vv'q)v) 



(vv, 



f\2 



(vv, 



>\2 



(5.24) 



Consider 



VVq _ , 1 ,, . 1- 



Vsi := {{x,v) E F+ I \x-Rv\ < R-S, \v ^v'q\ > 6, \v\ < -, d{x,v) < -}, (5.25) 
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— 1 1)1) 1 

Vs,i ■■= {ix,v)eF+\\x-Rv\<R-S--,\v--^v'o\>6+-, 



% 



\v\ 



<6-'-l-\ d(x,^)<^}, (5.26) 
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ioT < 6 < min(i?, Vq ^) and I e N, I > {R - 6y^ + 6. For < 5 < min(i?, Vg ^) and 
leN,l>{R- Sy^ + S, let X5,i e C^iF+) be such that 

Xv,^i<Xs,i<XVs,r (5.27) 

Finally, for < (5 < min(i?,V(^^) and m, / e N, / > (i? - 5)-^ + S, let 0^,^,/ G Co"^(F+) 
be defined by 

05,m/a;, t;) := x^/a;, t')p^ (t', t;) , (,_,,)(„/_(,„>),) • (5.28) 

(See (15:22|) . dSTM]) and flOTD - flOSD given below.) Note that from fl5:28|) and flOTj) it 
follows that 

SMY)Y>(t>S,m,l C Vs^i. (5.29) 

Using (EUl), (15:291) and (ESSl), it follows that 

lim /i(05,m,/,e) = 0, (5.30) 

for < 5 < min(i?, v^^^) and m,l eN, I > {R- Sy^ + 5 (we used that (xg + i^-Ug, Vq) ^ 

Note that using (15.251) we obtain 

XvsA^,v) < Xg{x,v)<1, leN, l>{R-6)-' + 6, (5.31) 

}^'^XVsi{x,v) = xvsix,v), (5.32) 



Z^oo 



for (x, v) G F+ and < (5 < niin(i?, Vq ^), where G is the compact subset of F+ given by 

G = {{x,v) e F+ \ \x - Rv\ < R - 6, \v - ^v'q\ > 6, \v\ < r^} and 



1)1) 

Vs := {ix,v) eF+\\x-Rv\<R-6, \v - ^v'q\ > 5, \v\ < 5'^ d(x,v) =0}. (5.33) 
Note also that, as n > 3, we obtain 

Xvs ix + Ri), v)dxdv = 0, < 6 < min(i?, v^^). (5.34) 




'V Jn^(R) 
From (13161) . (E29]), (15:3TD -( 1534D . it follows that 



limsuplimsup \l3{(f)s,m,h£)\ < s' for any e' > 0. (5.35) 
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Hence 



lim limsup \l3{<Ps,m,h£)\ = 0. (5.36) 



Let < 5 < min(i?,V(^^) and m e N. Using (^28\f . (ISlMj) and IKT3\f . we obtain 



l2i<Ps,m,l)= f f f5,mAt'^v)dt'dv, (5.37) 

Jv J-R 



for I > {R — 6) ^ + 6 where 



Pol 

X [xsA^ + -R^i w)£^(a;, t, V, v'q)] t=t(.'^y^,t',.) , (5.38) 

for t' G M and v ^ V such that i; and v^ are hnearly independent, where {t{x'(^,VQ, 
t', v), x{xq, v'q, t', v)) is defined by fl3.15p for v e F and t' e M. 

Using the estimates cr >0, \v\~^ < Vq^, < X5,i < 1, we obtain 

\f5,m,iit',v)\ < -{k + ~k){x'Q + t'viv'Q,v)), (5.39) 

Vo 

for I > [R — 6)^^ + d , f E M. and v E V such that v and v'q are hnearly independent. 
Using (12.91) and (13.11) . we obtain that the function arising on the right-hand side of 
fl5:39D is integrable on ^ x {-R, R). In addition from (I5:25|) - (l5:27l) . (15321) and (lOSj) . 
it follows that 

fs,m,ii^', v) -> fs,mit', v) as I ^ +00, (5.40) 

for t' eM. and v E V such that v and v'q are linearly independent, where 

Am(t',w) := ^(A;-^)(x[, + tX,^;[„^)p„(t',^;) (5.41) 

Pol 
X [XVs (x + Rv, v)E{x, t, V, v'q)] t=ti.'^,.>^,t' ,.) , 

x^x(^x^,,v^,t ,v) 

for t' G M and v eV such that v and v'q are linearly independent. Therefore we obtain 
by the Lebesgue dominated convergence theorem that 

lim ll{h,mA= I I f5,Ut',v)dt'dv. (5.42) 

'^+°° Jv J-R 

Let < 6 < min(i?, Vq""*^). We also have fs,m{t', v) < VQ^{k + A;)(xq + I'v'q, v'q, v)), for 
meN. From ([532]), (ICTD and (l5:22D . it follows that 

lim lim /a' (05,^,0 = / / ^1^ - ^K^^o + ^X, ^o,^) (5-43) 



m— »+oo /— ►+00 JyJ^fi \Vq 

[Xvsix + Rv, v)E{x, t, v, v'q)] ,^_(,, +,,„>), dt'dv. 

x — tv — xL-\-t'vL 
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From (15.431) . we deduce 



S—*0+ m—>+<X) l^+oo Jyj-fi |(^ol 




R 1 



lim lim Mm I2{(p5,ni,i) = ml^ ~ ^K^o + ^X^^O' ^) (5-44) 



X E{x, V, t, v'q) ,^_,,,+,,j., dt'dv. 






From (I3AD . it follows that 

9 f? ^ 

\ll{(t)&,m,i)\ < — ll^p(a;o + i''^0'Vo)IU°°(Kt/) sup {E- E){x,t,v,v'Q) , (5.45) 

Vq (a;,tj)eF 

tSM 

for < 5 < min(i?,V(7^), m G N and / G N, / > (i? - (5)"^ + 5. Using ([33)) and (IXTOll . 
we obtain 

|-^2(0<5,m,/,e)| < M--4||£(Li(F_),Li(i7'_^)) + \Ii{(j)s^rn,ue)\ + |/3(05,m,/, e) | , (5.46) 

for < (5 < min(/2, vc^^), m G N and / G N, / > (i? - 5)-^ + 5. From (15:461) . (EH]) and 
( 13A211 . it follows that 

|/2(^5,m,OI < M-^|U(Li(F_),Li{^+)) + l^2(<^'5.™,/)| + I lim /i((;/)5,„,/, s) | 

+ limsup|/3((^5,m,i,£)|, (5.47) 

for < (5 < min(i?, V(^^), m G N and / G N, / > (i? - 5)-^ + 5. Estimate fl319l) follows 
from ([52ZD, dSiD, (15361) . fl3:T2|) . fl5:4iD and (15^51) . D 



Proof of Theorem 3^. The method we use to prove (I3.22p is the same as in [13]. Let 



(a, A;), {a, k) E M.. Let f = a — a and consider Pf the X-ray transform oi f = a — a 
defined by P/(x, e) := /+^ /(t^ + a;)rft for a.e. (x, ^) G T§"-i. 
From dSI]) and /|x e H'^+^X), it follows that 



^-i(^)<A(n,X)||P/|U (5.48) 



where 



\\Pf\U--=( f I \Pf{xM''dxde\ 



and Di{n, X) is a real constant which does not depend on / and lie := {x G 1^" \ x9 = 0} 
for 9 G §"-^ Using ([SII]) (and (a, A;), (a, fc) G M), it follows that Pf{x,e) = for 
(x, 6*) G T§'""^ and |x| > /?. Therefore using also (15.481) we obtain 

^-i(^) < /^2(r^,X)||P/|U»(TS«-l), (5.49) 
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where D2{n, X) is a real constant which does not depend on a. 
We also use the following interpolation inequality: 



H'^iX) 



< 



2a + l 



H^^ 



n + l + 2f 



for -i < s < f + r. As a G M, it follows that 



(5.50) 



kiloo < D-i{n,f-)\\a\\^n^, < D^{n,f)M. 



Therefore, 



R 



a(xo - sv'Q)ds < 2RDs{n, r)M, 



for a.e. {x'o,v'o) E TS''~\ From fl532D it follows that 



-kol ^ I-R'^i^O-^'"0''"o)'^'^ — p~H\ Vifl'^K-^O'^'o)'^* 



-2V(7iiJD3(n,f)M 

> \Pia-a)ix'„v',)\, 



Vo 
for a.e. {xq,Vq) e M" x V^, Xqv'q = (we used the equality e' 



(5.51) 



(5.52) 



tl (,t2 



(5.53) 

e%t2 - h) for 

ti < t2 E M. and for some c E [ti, ^2] which depends on ti and ^2)- (In fact, the estimate 
( 15331) is valid for any (x^, v'^) eWxV, x'qv'q = 0, such that {x'q + tv'Q\tER}nXj^0 
or {x'o + t?;; I t G M} n X = 0.) Combining fl533|) . dEH, and (KTE^ . we obtain 



D2{n,X)V, ll^-^ll/z-^W ^ P-^IU(LHF-),LMF,)). 



(5.54) 



Combining (1534]) and (l530D . we obtain (13:221) . 

We now prove (13.231) . Using \vq\~^ > V^^ for v E V, (13. 2p . and (I5.5ip . we obtain 
that 




l^ol' JyJ-R 






> 



-4v~'^RD3(n,f)M 




I A; — k\{xQ + t'vQ, v'q, v)\dt'dv, 



(5.55) 



V J~R. 



for any (xq, v'q) G M" x V^, Xq^q = 0, such that x'q + sVq G X for some s G M, and where 
(^(xq, Wq, t', t;), x(xo, v'q, t', v)) is defined by (I3.15P for u G F and t' E {—R, R). 
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As {a, k) G 7VI we have ||o"p||oo < M. Using the latter estimate, (13.21) . and \v\ ""^ < Vq ""^ 
for all i; G y, we obtain 

||ap(x'o + tX,t;[,)|U^(M^,) sup \E - E\{x,v,t,v',) < Mc^^o'^d.^m 

ten 



X sup 

(x,v)SF 



1 /•* 1 pR+{x-tv)v'Q 

' a — a\[x — sv,v)ds + —j- I \a — d\{x — tv — sv'Q,v'Q)ds 



V\ .l-R 
< 4Rv-'Me'^^o'RD,in,r)A4. _ ~ 



V, 



01 ^0 



(5.56) 



for any (xq, iQ G W^xV, Xqv'q = 0, such that Xq+I'v'q G X for some t' G M. (We also used 

|e" — e"| < e™^^*^'"'''"'^|ii — n| where u = —\v\~^ J_j^a{x — sv,v)ds—\vQ\^^ f^ "*" ^ ° a{x — 

til — sv'q, VQ)ds and u denotes the real number obtained by replacing a by a on the right- 
hand side of the latter equality which defines u; using (I5.5ip (for a and for a) we obtain 
max(|n|, |{t|) < ARvQ^Ds{n,f)M.) Note that ||cr — cr||oo < Ds{n,r)\\a — d-\\^^+r for 
< r < f (see (I5.5ip ). Therefore, combining (I5.55p . (I5.56p . (I3.19P and (I3.22p . we obtain 

Let us finally prove l^2M- Let < r < f and let 9 = ^^'~_^^- . From <K2^ it follows 
that 



R 



RJV 



k{xo + tX, v'o, v) - Hx'o + t'v'o, Vq, v) 



dvdt'dxQ 



< D^WA - ^||£(Li(F__),Li(^+)) {} + 11-^ " •^ll£(Li(F_),Li(F+))j 



(5.57) 



where D4 = C2 Jx' v' =0 dxQ and C2 is the constant that appears on the right-hand side 

Wo\< R 

of (13.231) . From (13. ip . (I5.57P and the change of variables "x = Xq + t'v'f^" , it follows that 



"xV 



Hx,v'q,v) -k{x,v'Q,v] 



< Di\\A - '^\\c{L^F-),L^F+)) 



dvdx 



\A- A\\^-^ 



Integrating on Wq G V both sides of (I5.58p . we obtain (I3.24p . 



6 Decomposition of the albedo operator 

We now prove Lemmas 12.71 and 12.81 



(5.5^ 



D 
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Proof of Lemma 2.7. Let ip- G L^{0, \v\dxdv). Using the definition of K, we obtain 



(7^^_)(x,^;) = {K'ij^).^Jx,v) 



F+ 
1 1 



(6.1) 



2R />R+{x-tv)vi 



k{x — tv, Vi, v) 



VxV rl l"^!! Jo Jo 
xk{x — tv — tiVi, v' , Vi)Eq{x, v,x — tv,Vi,x — tv — tiVi) 

xiIj_[x — tv — tiVi, v')dtidtdv'dvi, 



for a.e. {x,v) G F+, where 



Eo{x, V,X-tV,Vi,X-tV- tiVi) = e-l''!"' ^O '^{x-sv,v)ds-\v,\-^ /„l a{x-tv-sv\,v,)ds ^ ^g^2) 

for X G M", t, ti G M and v, Vi G V . We recall that 7^ is a bounded operator from 
L^(P,\v\dxdv) to L'^{F+), i.e. 



||7^V^||F+ < C|||v|^||o, for any ijj G L^{0, \v\dxdv). 



(6.3) 



Hence we obtain, in particular, that the integral in t, ti, v' and wi, on the right-hand 
side of (16. ip is absolutely convergent for a.e. {x,v) G F+. 

Let us assume first that V = S"~^. Performing the changes of variables "x' = 
x — tv — tiVi" ("(ix' = t'^~^dtidvi"), we obtain 



(7^^_)(x,ti) = / /3{x,v,x ,v)tjj-{x ,v)dx dv 
'o 



(6.4) 



where 

/3(x, v, x',v') :- 



^^ ''k{x~tv,vi,v)k{x',v',vi) ^ ^ N 

; ; EniX, V,X — tV, Vi, X — tV — tlfl) 



(6.5) 



dt, 

t-^ ^\x — tv — x' I 

x — tv — x 
"1 = ET 



for a.e. (x,t;) G F+, (a;',!)') G O, where Eq is defined by (16. 2p . 

Now assume that V is an open subset of W^, which satisfies vq = inf^^y 1^;! > 0. 
From (16. ip . it follows that 



(7^V^_)(x,^;) = {K'ij_)^^Jx,v) 



hoo 
-1 / „n-2 



2R pR+{x-tv)uj 

r'^~^Xv'(^'^) / / k{x — tv,ru,v) 



X /i;(x — tv — tiu, v', ruj)EQ{x, v , x — tv , roo , x — tv — tiu) 
xiIj-{x — tv — tiuj, v')dtidtdrdv'duj, (6.6) 
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for {x,v) G F+, where for r > and a; G S" ^. Performing the changes of variables 
"x' = x-tv- tiw" ("rfx' = tl'^dtidu''), we obtain 



(7^^_)(x,^') = / p{x,v,x,v')'ilj-{x,v')dx'dv, (6.7) 



where 



f+oo r2R 

d(x,v,x',v') = 1— r / r 



1 r"" ._o /"^"^ . ,k{x-tv,ruj,v)k{x',v',ruj) 



n-2 
I 



Xv[ruJ) 



\x — x' — tv\"- ^ 
xEq{x, v,x — tv, ru, X — til — tiuj)]t^=\^-xi-ti\ dtdr, (6.8) 



ti LO^X — X' —tv 



for a.e. {x,v) G F+, (x',!)') G O, where i^o is defined by (16. 2p . 

From (E3D, (E3!)-(!63]), and (EZD-(EE1), it follows that for a.e. {x',v') G O, P{x,v, 
x',v') G L^(F+). Moreover from (16. 3p . it follows that the function O 3 {x',v') -^ 
P{x,v,x',v')i{j{x',v') G L^{F^) belongs to L^{0, \v\dxdv) for any i/j G L^{0, \v'\dx'dv'). 
Therefore 

\v'\-'l3eL°°{0,L^{F+)). (6.9) 

Now we prove ([232D- Assume A; G L°°(R" x T/ x 1/) and let 1 < p < 1 + ^, 
p'~ +p^^ = 1, be fixed for the rest of the proof of Lemma [ZTl We use (I6.10p . Using 
Holder estimate, the change of variables "y = x — tv" [dy = dxdt) and the spherical 
coordinates, we obtain 

2R ^ Y 

-^in-i ^^ ' dxdv 

Vs Jm+u^{R) \Jo \x — x' — tv\"- 






pfi)^// ^1^ < (2«)?Vo,(V,)Vo,(S"-.) (l^fl:-^, (6.0) 

JvsJ y^^ \y — x'y'^^ ^' n — [n — l)p 



\n—(n—l)p 
^ (y,v)eO 

for x' G M", Ix'l < 2R, where 

^5 := {v G 1/ I |i;| <r^}. (6.11) 

Assume first that V = S"~^ and let be a continuous function on F+. Then using 
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(I6.5p . cr > 0, Holder estimate and fl6.10p (with 6 = |), we obtain 




(x, v)P{x, V, x', v')dxdv 



V JRv+U^{R) 



<\\k\\l[ I I \<p{x,v)fdxdvY 



< \\k\ 




p{x,v)\ 

VJRv+UyiR) Jo 



^^ dtdxdv 



\x — tv — X' 



J\n-l 




X 




VJRv+n^{R) J 

2R 



V JRv+n^{R) \Jo 



\x — tv — x'|"~^ 



dt I dxdv 




<C{ I I \(P{x,v)fdxdv]' 

'v JRv+n^iR) J 



(6.12) 



where C= (2i?)^||A;||^Vol(§ 



n-l^ 



(4R) 



n — {n — l)p 



which proves 02.321) for V^ = S' 



n-l 



n—{n—l)p I ' 

Now assume that V is an open subset of M" which satisfies inf^gy \v\ > 0. Let e' > 
an 5 > be positive real numbers. Let be a compactly supported and continuous 
function on F+ such that supp0 C {{x,v) G F+ | \v\ < S~^}. We use the following 
lemma, whose proof is postponed to the end of this section. 

Lemma 6.1. The nonnegative measurable function (3i defined for a.e. {x,v,x',v') G 
F+xO by 



(3i{x,v,x',v') 



1 

\v\ 



-oo r2R 



Xv[ruJ^ 



k{x — tv, roo, v)k{x', v', ruj) 
\x — x' — tvl""^ 



X Eq{x, v,x — tv, ru, X — tv — tiu)] tj^=\^-x'-tv\ dtdr, (6.13) 



belongs to L°° {O , L^ (F^)) , where Eq is defined by (16. 2p . 
Let Ms' > vo be defined by 



Me' = Vo + £ 

From (16. 8p . it follows that 



/-I 




f]i{x,v, x' , v')\v\dxdv 



V JRv+U^{R) 



L°°{0) 




V JRv+U^iR) 



(6.14) 



{x, v)\v\f]{x, V, x', v')dxdv = Ii{x' , v') + h^x' , v'), (6.15) 
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for a.e. {x', v') G O and where 



h{x\v') 



(p{x,v) 



(6.16) 



Vs JRv+n^{R) 



VO 







Xv{ruj) 



k{x — tv, rcj, v)k{x', v', ruj) 



\x — x' — tv 



n-1 



Hx'.v' 



xEo{x, v,x — tv,ru,x — tv — tiu;)]^ ^=x_x'-tv dtdrdxdv, 

(j){x,v) (6.17) 

k{x — tv, ruJ, v)k{x', v', roj) 



Vs JRv+n^{R) 



hoo /.2_R r 

„n-2 



M, 



XvVruj) 



\x — x' — tv 



n-1 



xEq{x,v, X — tv, ruj, X — tv — tlUJ)]^ ^=x_xi-tv dtdrdxdv. 

Using fl6.17p and the estimates r"~^ = r^-'^r~^ < Mj^r""^ ioi v E V and r > M^/, and 
using (16.141) . we obtain 



veV JRv+U^(R) 

< £:'||0||loo(f+), for a.e. {x',v') G R 
From (16.161) and Holder estimate, it foUows that 



Pi{x, v, x" , v")\v\dxdv 



L°^{0) 



(6.18) 



\h{x',v')\ < 



Vs JRv+U^{R) Jvo 



M^i r2R 

r"~^dr 







\x — x' — tv 



I — +oAn-l 



-dtdxdv. 



2R 



1 



< II A; II ^^^^^' ^ ^^_^. 

n — 1 \.lvs J Rv+u^iR) \Jo \x — x' — tvr 



dt I dxdv 



X I / / \(f){x,v)\'^ dxdv 

'vs JRv+n^iR) 



(6.19) 



for a.e. {x' ,v') G O. Combining fl6J[5D . fl638D - fl639D and fl630D . we obtain fl232D with 



C{e',5,p) = \\k 



,^^ -i2R)7 {Vo\iVs)Vo\iS-~')y [^ 



\2 ^'-"e' 



n 



n — {n — l)p 



D 



Proof of Lemma \2. <g[ Let 0_ G Cq(F_) (which denotes the spaces of C^ compactly 
supported functions on F_). Let := (/ + K)~^J(f>_. Then note that 



J(/)_ - KJcf). + K\I + Ky^Jcj). 



(6.20) 
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Thus 

A<P. = (j)\F^ := {J(t>-)\F^ - {KJ(t)-)iF^ + 7^(/ + Ky'j<P_. (6.21) 

Prom (Km and ([OSD, it follows that 

(J0_)ip (x,v) = / / ai(x,v,x',v')(l)-(x' — Rv',v')dx'dv' , (x,v) G F+. (6.22) 

•'^ •' W\<R 

From the definitions of K and J, we obtain 

|a;'|<-R 




{KJ(f)^),p {x,v) = / / a2{x,v,x',v')(j)_{x' — Rv' ,v')dx'dv', {x,v) ^ F^. (6.23) 



Lemma EJ] follows from (lOB - ( 16:231 ) and Lemma O □ 

Proof of Lemma \6.1\ Using (16.131) and the estimate a > and using the change of 
variables "y = x — t{)" {dy = dtdx) and (12.91) . and spherical coordinates, we obtain 

/ / \v\(3i{x,v,x',v')dxdv 

+00 2R 



< I I I r"^i 



Xv[ruJ^ 



k{x — tv, ruj, v)k{x', v', ruj) 



\x — x' — tv 



n.-l 



dtdrdxdv, 

x — tv — x' 



V Rv+Uv{R) vo \x-tv-x I 

+ 00 

-r-, [;^ / r"-'^ [Xv{ruj)k{y,ruj,v)k{x',v' ,ruj)]^^ y-x' drdydv 

\X y\ J \y-x'\ 

V ySK" vo 

\y\<R 

+ 00 

X\y\<Ri^' + 1^'^) / ^^ ^Xv('"^)c"p(2;' + r'u, ruj)k(x\ v', ru)drdr'diO 

J 

S"-i vo 

+ 00 

< ||o-p||oo-R / / r'^~^xv{ru!)k{x',v',ru)drdu; = R\\ap\\ooCrp{x',v') < R\\ap\\l^, 

§"-1 vo 

for a.e. {x',v') G O. The lemma is proved. D 
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7 Proof of existence of the albedo operator 

In this section, we prove Lemma 12.41 and Propositions 12.51 and 12.61 

Proof of Lemma \2.4\ Using the definition of Ti~^, the estimate a > and (12. lip , we 



have 

IlklTrVllo < / / / / \f{y+{w-t)v,v)\dtdwdydv 

Jv Jn^{R) J-R Jo 

p p pR pw 

= / / \f{y + tv,v)\dtdwdydv<2R\\f\\o, 

Jv JUviR) J~rJ~r 

for / G L\0). 

Using the definition of A2 and (12.111) , we have 

\\A2\v\-'f\\o < f [ [ k{y + tv,v',v)\v'\~'\fiy + tv,v')\dtdydvdv' 

JvxV JUuiR) J-R 

= 11 [ My' + t'v',v')\v'nf{y' + t'v',v')\dt'dy'dv' 

Jv Jn^i(R) J~R 

< llk'r^^p(a;',^^')IU-(o)||/||o, 
for / e L^{0). We also used (12.91) and the change of variables 

p /'+OO p p-\-oc 

/ / f{y + tv)dtdy= / f{y' + t'v')dt'dy\ (7.1) 

J yV^O J —00 Jy'v'^0 J —00 

for /€ L^(M") ajidv,v' e V. 

Using the definition of Ti~^ and A2 and Lemma [2. 11 (I2.23P and (17.10 . we obtain 

p 1 /" r^ pR+{y+wv)v' 

< —- / k{y + WV,v',v) / ^-\v'\'^f^^a^(y+wv-sv',v')ds 

Jvxv P'l Jn^(R) J-R Jo 

X \f[y + wi) — tv', v')\dtdwdydvdv' 

-| p pR pR+w' 

±. / k{y' + w'i',v',v) g-K|-^/>p(?/'+K-^Ky)rf^ 

VxV l"^'! Jn^fiR) J-R Jo 

X \f{y' + (w' - t)v', v')\dtdw'dy'dvdv' 
{--^e-\'''^'"^^^'''-'-y'+''^'''''^'''dw']\fiy' + ti/,v')\dtdy'dv' 




'V JU^,{R) J-R \Jt dw 
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1 _ g-|.rVo''-^"'-p(-+-'y)d. j \f(x',v')\dx'dv' < (1 - e"2^^o-'lkpll-)||_^ no, 

for / G L\0). 

Item iii follows from items i and ii (under (12.241) . we also use that ||y42Ti^^|| < 

Proof of Proposition l27R We first prove item i. 
Assume ( lOTj) . For all / G D(T), 

T/ = (Ti + A2)/ = (/ + A2Ti-i)Ti/. (7.2) 

From (I2.27P it follows that T admits a bounded inverse in L^{0) given by T~^ : = 
Ti~^(/ + y42Ti~^)~^. Using the latter equality, we obtain 

(/ + K){I- T-'A2)=I + T^~'A2 - Tr^(/ + A2T^^yA2 

-Ti-\I + AaXr^ -/)(/ + A2Ti-yA2 = I. (7.3) 

The proof that (/ — T'^A2){I + K) = / is similar. We now prove that (12.261) implies 
(II2ZD. For / e D(T), 

T/ = (Ti + A2)/ = Ti(J + T1-M2)/ = Ti(/ + JO/. (7.4) 

Let us prove {I + K){D{T)) = D{T). From the latter equality and (17^ it follows that 
T admits a bounded inverse in L^{0) given by 

T"^ = (/ + A')~^Tr^ (7.5) 

As K = Tr^A2, we have (J + K){D{T)) C D{T). Let c/ e L'(T), and let / = 
(/ + K)-^T^-'g G L\0). Then / = -KT^~'g + g = -T^~'A2T^-'g + gE D{T) (we 
recall that g G D{T)). 

Equality (O) still holds. Using (ED, (ESD and the fact that Ti : D(T) -^ L^{0) 
is one-to-one and onto L^{0), we obtain (I2.27p . Item i is thus proved. Item ii follows 
from item iii of Lemma [2.41 and item i. We shall prove item iii. Note that (see (17. 3p ) 

(/ + K){I- T-M2) = / = (/- T-'A2){I + K) in C{L\0)). (7.6) 

Note also that L^(0, \v\dxdv) C L^{0) and recall that K is a bounded operator in 
L^{0, \v\dxdv). Therefore, we only have to prove that T~^y42 defines a bounded operator 
in L^{0, \v\dxdv). Note that 

T~^ = Ti-\I + A2Ti-^y\ (7.7) 

From item i, (17. 7p and item i of Lemma [2.41 it follows that T"^A2 defines a bounded 
operator in L^{0, \v\dxdv). Thus item iii is proved. D 
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Proof of ProposiUonlKE Let /_ G L\F_). From (12191) . it follows that J/_ G W. 
Hence J/_ G -^^(O, |f |(ix(it;) and from (12.280 it follows that ( 12.25^ is uniquely solvable 
in L^{0, \v\dxdv) and its solution is given by (/ + K)^^J f^ which satisfies 

\\{I + Kr\jf4\L^oMd.dv) < C^o||/-||f_, (7.8) 

where Co = ^R^i + v'^^\\a\\oo)\\{I + Ky^'W^moHdxdv))- 
Let f := {I + K)~^Jf^. Hence by definition 

f = Jf^-Kf mL^{0,\v\dxdv). (7.9) 

Using (17.91) . we check that the following equality is valid in the sense of distributions: 

To/ = -A^f - A^f. (7.10) 

Using (moll , we obtain To/ G L\0) and 

\\Tof\\ < (Iklvr^lloo + \\\v\~'^(rp\\oo)\\f\\LHO,\v\dxdv)- (7.11) 

Therefore f E W (item i is thus proved), and using ( 12.141) and ( 17. lip we obtain 

||/|f+||f+ < max((2i?)~\l)(|||i;|~V||oo + |||t^rVp||oo + l)\\f\\Li{o,\v\-^dxdv)- (7.12) 
Item ii follow from (EHD and (HA^ . 
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